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The "recoil" correction of order ma to the hyperfine sphtting of positronium 
ground state was found. The formahsm employed is based on the noncovariant per- 
turbation theory in QED. Equation for two-particle component of full (many-body) 
, wave function is used, in which effective Hamiltonian depends on the energy of a 

\^ I system. The effective Hamiltonian is not restricted to the nonrelativistic region, so 

' there is no need in any regularization. To evaluate integrals over loop momenta, they 

. are divided into "hard" and "soft" parts, coming from large and small momenta re- 

Q I spectively. Soft contributions were found analytically, and hard ones are evaluated by 

C3 ' numerical integration. Some soft terms due to the retardation cancel each other. To 

calculate the "hard" contributions, a great number of noncovariant graphs is replaced 
Qh, by only a few covariant ones. The hard contribution was found in two ways. The first 

way is to evaluate contributions of separate graphs, using the Coulomb gauge. The 
second one is to calculate full hard contribution as a whole using the Feynman gauge. 
The final result for the "recoil" correction is 0.381 (6)ma^ and agrees with those of 
||l|, ^. Diagram-to-diagram comparison with the revised results of Q was done. All 
^ ■ the results agree, so the "recoil" correction is now firmly established. This means a 

. considerable disagreement with the experimental data. 



1. INTRODUCTION. 



The treatment of relativistic bound states is one of challenging tasks of Quantum Elec- 
trodinamics. Positronium, the bound state of the electron and the positron, is one of the 
most appropriate objects for theoretical and experimental study of relativistic bound states. 
Because of the small masses of its constituents, effects of the strong and the weak interac- 
tion are negligible compared with the accuracy of current experiments on the positronium 
spectroscopy. On the other hand this experimental accuracy is good enough for compari- 
son with results of modern theoretical investigations, which reached accuracy level ~ ma^ 
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for contributions to the energy. Therefore, now positronium spectrum can be investigated 
within QED framework. 

Now the most accurately measured positronium property is hyperfine sphtting of the 
ground state, i.e. energy difference between the l^Si and l^S'o-states, denoted below as 
Ai/. The two best experimental results for this value are 

Au = 203 387.5(1.6) MHz , (1) 

Au = 203 389.10(0.74) MHz ; (2) 
obtained in |[ and |0 respectively. 

Calculation of Au has a long history; contributions ~ ma^, ~ ma^, ~ marina were 



found in [|^ H, ^ [10], |Tl|, [T^ ; their sum equals 

^ ^ ^ = 203 m.29MHz . (3) 
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To compare experimental results (|],0) with the theoretical one, the contributions ~ ma^ 
(without the logarithm) are to be obtained. An essential progress in their calculation was 
made only recently. 

There are several sets of corrections ~ ma^ to the hyperfine splitting, that are different 
in origin and may be calculated independently. First, these are three sets consisting of 
contributions arising from the one-, two-, and three-photon annihilation; they were found 
T3| , |1^], ||15| and [|l^, |l^ respectively. Second, there are contributions associated with 



m 



radiative corrections to the Breit potential, i.e. those of formal order (Za) a m), found 
in 



[T^. Third, there exist "recoil" corrections, i.e. those coming from the graphs in 
which all photon lines connect two fermion ones; these corrections are of conventional order 
(Za)" m (here Ze is the charge of one of the particles; Z = 1 in positronium). Finally, 



there are radiative-recoil (i.e. ~ (Za) C(m) corrections found in 20, 21, 22 



Below the sum of contributions ~ m{Za)^, ~ m{Za)^ In a to Au is denoted as Az/^ec- 
Calculation of its nonlogarithmic part is the most difficult step of calculation of Au. This 
was the subject of works [l|, ||, |^. The present work is also devoted to calculation of 

The results of the three first works, namely |T], ^, were all different. The result 
of work coincides essentially with the result of [0] (they are equal to 0.3763ma^ and 
0.3767(17)ma® respectively). The result of the present work is 0.381(6)ma® and also coin- 
cides, within its accuracy, with the results of [jl], 0. Recently the result of was corrected 
by its authors. Now it is equal to 0.3764(35)ma^ and agrees with the results mentioned 
above. 

In the so-called NRQED (nonrelativistic QED), an effective field theory equivalent 
to QED, was formulated and then applied to calculation of Aurec', nonlogarithmic part of 
the result equals 0.167(33)mQ;^. However, just now the preliminary result was obtained in 
the NRQED framework |^^, which is consistent with that of [|^, 0. In the calculation 
was performed by Bethe-Salpeter formalism. In |jl]] the calculation was done by an effective 



2 



Hamiltonian approach; the method employed in |0] in fact also uses an effective Hamilto- 
nian. Difference between |[| and 0] is in different regularizations used; besides that, the 
result of HI] was obtained by numerical integration, whereas that of is analytic. The 
effective Hamiltonian approach is essentially a combination of ideas of the NRQED and 
the old-fashioned noncovariant perturbation theory in QED. The latter was used first for 
QED bound state calculations in ||2^ . 



The present paper also employes a formalism based on old-fashioned noncovariant 
("time-ordered") perturbation theory for QED. Starting from the Shrodinger equation for 
full (many-body) wave function one easily obtains an equation containing only two-particle 
part of the wave function, in which effective Hamiltonian depends on the energy of a system. 
Then the effective Hamiltonian is divided into unperturbed part and perturbation, zeroth 
approximation for the two-particle wave function is found, and using them corrections 
to the energy levels are calculated by means of usual Rayleigh-Shrodinger perturbation 
theory. This calculation consist essentially in finding of expectation values of operators 
corresponding to graphs of the noncovariant perturbation theory over the "unperturbed" 
wave function. To evaluate the integrals over loop momenta each of the integrands is di- 
vided into two parts, "hard" and "soft"; here terms "hard" and "soft" mean that in the 
desired order a "hard" part entirely arises from the region of momenta ~ m, whereas for a 
"soft" part the region of momenta ~ am is also essential. This decomposition is performed 
in such a way that "soft" contributions may be easily found analytically, and "hard" ones 
are found by numerical integration. To calculate the "hard" contributions with the required 
accuracy external legs of the graphs may be set on the mass shell, which allows to replace 
the sum of a great number of noncovariant graphs by the sum of only a few covariant ones. 

The method of calculation described above differs from that of and in the way 
effective Hamiltonian is defined. In [|I| and 0, effective Hamiltonian is constructed so as 
to reproduce the scattering amplitudes, whereas in the present paper it is immediately 
derived from the full QED Hamiltonian. Besides that, the effective Hamiltonian used in 
the present work is not restricted to the nonrelativistic region, in difference with and 
0, so there is no need in regularization of the effective Hamiltonian or matrix elements it 
enters. 

The plan of the rest of the paper is following. In section 2, I briefly describe the formal- 
ism used in the paper, i.e. the effective "Shrodinger equation" for two-particle component 
of the wave function, "zeroth" approximation to the wave function, perturbation theory 
formulas for the calculation of Az/; these formulas are transformed into form of expecta- 
tion values over the nonrelativistic Coulomb wave function. In section 3, I describe the 
general method of calculation of these expectation values, namely, the explicit manner to 
divide them into "soft" and "hard" parts; correspondence between these values and co- 
variant graphs (and also scattering amplitudes on mass shell) is established. In section 4 
this general method is applied to contributions of various sets of graphs in turn. Section 5 
contains the description of the procedure of numerical evaluation of "hard" contributions, 
the checking of this procedure, and comparison of the results obtained with the results of 

0, 1^ . Section 6 consist of summary including comparison with experimental results. 
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2. FORMALISM. 



The starting point of the formahsm used in the present paper is the Shrodinger equation 

H\i;r^ >= E\^^ > , (4) 

where the positronium wave function \ipN > includes components with different number 
of particles. In (^) H is the full QED Hamiltonian, H = Hq + Vn , where Hq is the free 
Hamiltonian, is the interaction. 

Direct use of the many-body wave function is unconvenient. However in the lowest 
order in a the only nonvanishing component of > two-particle one, which is denoted 
as |e >. It is easy to obtain an equation containing only |e >. Let be projector onto the 
two-particle subspace, Pa = I — Pe, 

\e> = Peli'N > , \a> = Pal^pN > ■ 

From the Shrodinger equation one gets 

PeHPe\e > + PeHPa\a > = E\e> , (5) 

PaHPe\e > + PaHPa\a > = E\a> . (6) 
Expressing \a > from @ and substituting it into (|^) results in 

( P,HP, + Va{E) - E ) |e > = , (7) 



where 



Energy levels which could be obtained from ([7p are derived in a form of perturbative expan- 
sion. For this purpose "Hamiltonian" PeHPf, + Va{E) is divided into the unperturbed 
part = PeHoPe + Vo and the perturbation V{E) = Va{E) + PeVNPe - Vo , where 

Vo is the usual nonrelativistic Coulomb interaction. Contributions to the energy, which are 
of 1-st, 2-nd, and 3-rd order in V, have leading orders (besides logarithms) ma'^, ma^ and 
ma^ respectively. 

Clearly V{E) + Vq = Va{E) + PeVNPe is the two-particle irreducible kernel for the 
Green function appearing in the noncovariant perturbation theory (for transitions between 
two-particle states). Its expansion in a is 

ViE) + Vo = PeVN ( 1 + ^^Vn + Vn Vn + ) Pe 

and corresponds to the set of two-particle irreducible graphs of the noncovariant pertur- 
bation theory (here and below "irreducible" is meant in noncovariant sense, for instance, 
graph shown in fig.l is irreducible). 

Since diagram technique for the noncovariant theory is not of common use, its rules are 
described here. External fermion lines correspond to positive-energy spinors normalized as 
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W^u = 1 (for convenience, the fermion-fermion channel, instead of the fermion-antifermion 
one, is considered). Internal fermion lines are described by the projectors 



A±(p) = [Ep ± (70 + ap))/ (2£p) , 6p = ^p2 + m2 
factors corresponding to external lines may be written in form 

uip) = ((25,)/(e, + m))V2 A+(p) 




is (2 X 4) matrix, w is a two-spinor normalized as w~^w = 1 . The two fermions 



are assumed to have charges of opposite sign; then Coulomb photon with the momen- 
tum q gives the factor —Ana/q^^ , and magnetic (transverse) photon gives the factor 
( —Ana ) ( did2 — (aiq) (tt2q) /q^ ) /2g , (in the noncovariant perturbation theory en- 
ergy of virtual photon is equal to its momentum). Every negative-energy projector A~ 
gives extra factor (—1). Factor { E — Ek + iO )~^ corresponds to each intermediate state, 
where Ek is energy of the intermediate state, E is full energy of a system. Loop momenta 
qj which remain undetermined lead to the integration d^qi / (27r)^ . 

In the rest of the paper all operators and wave functions refer to the two-particle sub- 
space. All the calculations are performed in the center-of-mass frame, so the only variable 
except spins is p, i.e. momentum of particle 1 in the c.m.f.. Eigenstates of p are normalized 
as usual, according to < p'|p > = {27t)^ <5^(p ~ p') , kernel of an arbitrary operator X is 
denoted as X(p, p') = < p'|X|p >. 

The kernel ^(p, p'; E) + Vo(p, p') is given by the sum of matrix elements corre- 
sponding to all the irreducible graphs for initial and final states having momenta p and p' 
respectively. The tree diagram with Coulomb photon corresponds to operator V^, and that 
with magnetic photon corresponds to Vm- Let us write V in the form 

V{E) = Vi + V2 + Vs + , (8) 

where 

Vl = Vi, + Vm , Vi, = K - , 

V2, V3 and so on correspond to sums of the irreducible graphs having two, three, and more, 
photon lines respectively. Being expanded in f/c up to the second order, H^^^ + Vi gives 
the Breit Hamiltonian. 

Consider, for example, the first-order in V correction to energy levels. The first term 
in (P) gives to this correction contribution of formal order ~ ma^ , the second term gives 
~ ma^ , and so on. However in fact this expansion does not converge: for instance, spin- 
independent correction to the energy, arising from V2, V3 and further terms of expansion 
(I), are all of order ~ ma^ . This takes place due to graphs shown in fig. 2: (a), (b), (c), 
and so on, because these graphs produce "ultrasoft" contributions, i.e. those coming from 
the region where momentum of the virtual magnetic photon is ~ mo^ . Fortunately for 
the hyperfine splitting the "rest term" of this expansion corresponds to the graphs first of 
which is graph shown at fig 2(c), and has order ~ ma^ . 
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When deriving formulas for corrections to energy levels the ground state of positronium 
may be regarded as non-degenerate, because S-states with cr = and a = 1 (where a is 
the positronium full spin) do not mix. For perturbation V depending on full energy of a 
system, correction to non-degenerate energy level number n reads 

IT/ 12 T/ T/ T/ 

A j-i Tr \ ^ I ^ nm\ \ ^ \ ^ ^ nm ^ mk ^ kn 

" = Vnn + 2^ ^ _^ + 2^ 2^ r E - E ) (E^ - E ) 
. y dVnn , ^ ( y Km? ] , ^ ( y2 9}^] .q^ 

where summation over continuous part of the spectrum, as well as over discrete one, is 
implicit; Em, En, Ek are zeroth-order energy levels; all matrix elements and their deriva- 
tives are taken at E = En. It may be easily shown that for calculation of Au to order ma^ 
it is enough to use only three first terms in (P). Keeping only terms of the order desired, 
the recoil contribution to Au equals (up to the order ma^) 

< ^1 ViG'ViG'Vi + + ViG'V2 + V3 + ViG'Vi + V2 + Vi\ij> |^=J , (10) 

where > is zeroth approximation to |e > , i.e. solution of equation 

i/W 1^ > = E2\iJ> , (11) 

and G' = G'{E2) is Green function of equation (|ll|) with the ground state pole subtracted 
out (the Green function without the subtraction is denoted as G). 

As H'^'^'> does not depend on fermion spins, > may be assumed to have the form 



(^V'(P) |P> \X> 



where |x > is spin part of the wave function. So to find > is essentially to find i 
i.e. to solve equation 

{T + - E2) = (12) 

where T is kinetic energy of the two particles (including the mass). 
Nonrelativistic approximation for (|12]) is 



in + Vo - Eo) Mp) = , (13) 
where Tq is nonrelativistic approximation for T. The ground state solution for (0) is 

Eo = 2m - ^ , (^o(p) = 8 j'/'n'/' ^ 
m 

(here and below notations 7 = am/ 2 , = + 7^ are used, for an arbitrary momen- 
tum k). 
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Let us write as 

V? = Co v?o + Ay? 



where Ay9 is orthogonal to (yJo! then ( p!2|) is rewritten as 

( To + K) - ^0 ) = ( - AT + ) , (14) 

where AE = E2 — Eq , AT = T — Tq . Solution of ([l^) may be written in the form 

= G'oiEo) {AT - AE)^ , (15) 

where G'q is Green function of equation (|TB[) (i.e. usual nonrelativistic Coulomb Green 
function) with the ground state pole subtracted out (corresponding Green function without 
the subtraction is denoted as Go)- 

Relativistic free two-particle Green function and nonrelativistic free one are referred to 
as 5*0 and S: 

So{E') = (E'-To)-' , 
S{E') = {E'-ry ; 
G and Go are expressed as expansions 

Go(^') = {E' - (To + Vo))-' = So + SoVoSo + , 

G{E') = {E' - {T + Vo)y' = S + SVoS + (16) 

Let us also denote 

T = Gg — 5*0 — S'oVqS'o . 
It may be shown by iteration of (|T5|) that with the sufficient accuracy its solution is 

^ ^ {1 + {S- So)Vo + SVoiS - So)Vo + LAT ) \e>=e, V>o ■ (17) 

It is convenient to use in the rest of the paper notations < > and < >(«) , defined so 
that for any operator X 



<X>^j y.+ (p')X(p,p',E)|^zJv,o(p) 



(27r)3 (27r)3 ' 

and < X >{n) stands for the sum of contributions ~ a", a" In a to < X >. 
If one rewrites (0) using 

G' ~ (S* + SVqS + L )\e'=Eo , 

expression for ip from ipo ([T7|), and keeping only terms of the relevant order, one obtains 

AUrec =< V3 + V2 + Uv2 + Uc + Um + UmM + UmCM + UmMM + Ul >(6) , (18) 

where 

Uv2 = {V2SV^ + V2iSV,-SoVo)) + h.c. , (19) 



and 



Uc = V^, + Uc2 + f/c3 , (20) 

Um = + Um2 + Umz , (21) 

Umm = Umm2 + Umms (22) 

Umcm = VmSVcSVm , (23) 

UmMM = Kn'S'KnS'Kn , (24) 

f/i = \/iL\/i + ( ViLAT + h.c.) , (25) 

Uc2 = Vi.SVi, + ( Vic{S - 5o)Fo + h.c. ) , (26) 

Um2 = Vm{SV, - SoVo) + h.c. , (27) 

UmM2 = VmSVm , (28) 

Uc3 = Vo{S - So)VUS - So)Vo + {Vi,SVo{S~So)Vo + h.c. ) + (29) 
+ ( Vi,SVi,iS - So)Vo + h.c. ) + VicSVoSVi, + Vi^SVi^SVi, , 

Um3 = {VcS - VoSo) Kn {SV, - SoVo) + ( Kn^K {SV, - SoVo) + h.c. ) , (30) 

Umm3 = ( K^5Kn {SV, - SoVo) + h.c. ) . (31) 



3. THE METHOD OF THE CALCULATION. 

In this section metliod of calculation of contributions to Ai/rec, written out in (|1^) - 
(0), is presented. The main idea of this method is just the idea used in other modern 
works concerning bound states, though there is some difference in its implementation. 
Every contribution to Az/^ec is divided essentially into two parts, "soft" and "hard" ones. 
Here the following is undermined. Contributions to Az/^ec are in fact integrals over loop 
momenta. Any contribution under the consideration is called "hard" if it is determined 
(with accuracy required) by region where all the loop momenta are of order m. Otherwise 
this contribution is called "soft". To separate "soft" contributions from "hard" ones, the 
integrands are to be expanded in powers of the momenta. In order to know whether given 
term of such expansion contain "soft" contribution, all the momenta should be set to 
~ am, and then simplest power counting let us know whether this momentum area give 
rise contribution ~ ma^ to Az/^ec (for the term under the consideration). If this is the case, 
the term considered clearly must be treated as "soft" contribution. On the other hand it 
can be proved that in the present problem absence of the contribution ~ ma^ due to the 
momenta region mentioned above lead to absence of "soft" contributions of the relevant 
order at all, at least for the way of calculation described in the present paper. 

It is clear that there is some freedom in explicit way of the decomposition into "soft" 
and "hard" parts, and "soft" part can always be defined so that its contribution to the 
integral is easily evaluated analytically. On the other hand, calculating "hard" part one 
can put Eq = 2m in the integrands, and the integrals obtained can be easily evaluated by 
means of numerical integration. 
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Contributions to Ai/^ec are naturally divided into three sets, different in form and 
calculation procedure used. These are "tree" contributions (< Vic, Kn >(6))! "one-loop" 
(< V2, Uc2, Um2, Umm2 >(6)), and "two-loop" ones (< V3, Uv2, Ucs, Ums, Umms, Umcm, Ummm 
(here words "one-" and "two-loop" imply that the operators inside < >(6) are integrals over 
one and two momenta respectively), and also term < Ul >(6) that corresponds to graphs 
having three loops or more; the last is determined by region of momenta ~ ma and can be 
easily found analytically. 

Consider different kinds of the contributions one by one. 

Each two-loop contribution has the form < X2 >(6) where 

^2(P,P',7) = / >"2(p,p',qi,q2,7) • (32) 

It is convenient to define the value n, the "divergency power" of the X2; namely, Y2 ~ 
at p ~ p' ~ gi ~ g2 ~ 7 ~ "C m . All the "two- loop" contributions have n > ; if n = 
then region gi ~ ^2 ~ contribute to order ~ ma^ and < X2 >(6) niay contain terms 
~ marina. It can be shown that for all the "two-loop" terms having n > the only 
essential contribution to the integral arises at p ~ ~ 7 , gi ~ g2 ~ |qi — q2| ~ ; de- 
pendence of Y2 on p, p' in this area can be neglected that results in 

< X2 >(6)~ |v5o(0)pX2|p=p/=^=o|^=o egw^f < X2 >p , (33) 

where 

MO) = j j^.Mp) =\/-. 

The integral in (|33|) does not depend on the small parameter a and so can be easily com- 
puted numerically. For the sake of brevity notation < >p defined in ( P5| ) is used throughout 
the paper. 

If n = the integrand I2 is to be divided into two parts: the soft part, denoted as 
Y21, and the hard one, F225 ^2 = ^21 + ^22, so that Y21, I22 have "divergency power" 
n = and n > respectively. It may be done in such a way that 1^21 equals approximately 
a uniform function at p, p',gi,g2,7 -C m (the only exception is the contribution of graph 
shown at fig.2(b), treated in the section 4). Similarly X2 = X21 + X22 where X21, X22 are 
integrals of form ( p^ ) with the Y2 replaced by 1^21, ^22 respectively. It may be shown that 
with such choice of 1^21 in all cases met in calculations below the correct recipe for the 
calculation is 

< X2 >(6) = < X2 - X21 >p + < X21 >(6) . (34) 

The decomposition into X21 and X22 is evidently non-unique, and X21 may be choosen so 
as readily to calculate < X21 >(6) analytically. 

In the method described above the < X21 >(6) is defined by the integral which is 
assumed to converge. Hence the 121 falls at gi, g2 m faster then uniform function having 
n = 0. As the scale of the momenta, at which the original function Y2 begins to fall rapidly, 
is gi ~ g2 ~ m, it is natural to demand the same behavior from the function 1^21 as well; in 
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other words, there exist effective "cutoff' of the I21 at gi ~ ^2 ~ Tn- This "cutoff' reminds 
regularizations in nonrelativistic effective theories; however, in the approach described here 
there is no need in artificial regulators such as large regulating mass or dimension different 
from 4. Similarly to what explained above, in the rest of the paper all functions looking 
like uniform (with n = 0) at small momenta fall rapidly enough at large ones. 

For calculation of "tree" and "one-loop" contributions it is necessary to find corrections 
to "leading" terms, which are of first and second relative order in the a. The method of find- 
ing of these corrections is rather obvious and is shown by the following simple illustration. 
Consider calculation of < X > where X = X{p, p') is sufficiently smooth function that 
doesn't depend on a, and expansion of X(p, p') — X(0, 0) in p, p' starts from third-order 
terms. Evidently 

<X>= {Ao + + A2)\:zl , 
Ao = |^o(0)|'X(0,0) , 

Ai = MO) I ^</^o(p) (X(p,0)+X(0,p)-2X(0,0)) , (35) 

A2 = J (^^0 (pO (^(p, p') - ^(P, 0) - X(0, p') + X(0, 0)) ^o(p) . (36) 

It is easily seen that when calculating three first terms of expansion of < X > in a, in (^,^) 
V9o(p) may be replaced by ; hence ^0,1, 2 exactly have orders , , 

respectively. 

Each "one-loop" contribution has the form < Xi >(6) where 

^i(P,P',7) = / ^a'n(p,p',q,7) ; (37) 

Let n to be defined so that Yi ~ 5"""^ at p ~ p' ~ g ~ 7 ~ 5 <^ m (if one writes Xi^qVo 
as the loop integral, the n is its " divergency power" ; n = imply that the region q ~ ma 
contributes to order ~ ma^, and may cause terms ~ marina) to exist. 

Let us write Yi for any given contribution as sum of "soft" and "hard" parts, Yi = Yn + Y12 , 
where Yn and Y12 has n < and n > respectively (similarly Xi = Xn + X12). It may 
be done so that Yu at p,p',q,'-f <^ m equals approximately a uniform function (the only 
exception is the contribution of graph shown at fig.2(a), which is treated in the section 
4). It can be proved that under such choice of Yu contributions coming from Yu can be 
obtained, for all terms to treat, using formula like (p5D, which results in 



<Xi>(6) = <((Xi-Xn) - (Xi-Xn)o) W + h.c. >p + < Xn >(6) ; (38) 

in ( p8D the notation ( )o , is used, defined so that (X)o(p,p') = X(p, p')|p=p/=o • Evi- 
dently Yii can always be choosen so that second term in (pH]) is easily evaluated analytically. 



"Tree" contributions to Aurec are calculated in the way just like one described above. 
There are only two "tree" contributions to be found: < Vm >(6) and < Vic >(6)- Each of 
them is divided into "soft" and "hard" parts (the explicit way of this partition is described 
in the section 4), and then the "hard" contributions are evaluated according to (^61). 
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In the section 4 separation of "soft" contributions from Ai/^ec is done using formulas 
that results in 

< >(6) = < 1^3 - Wv2 >p + , (39) 

<V2 + Uv2 >(6) = < V2SVm + h.c. >p + < (^2^14 - Wv2) + h.c. >p + E^^ , (40) 

< Uc >(6) = < V,SV,SV, -Wc>p + Eg , (41) 

< Um >(6) = < V.SV^SV, - Wmi >p + < {Vrr^SV.SV, - Wm2) + h.c. >p + El, , (42) 

< Umm >(6) = < {VmSVmSVc - Wmm) + h.c. >p + Efj^ , (43) 

< Umcm >(6) = < ymSVcSVm — Wmcm >p + E^Wj^ , (44) 

< UmMM >(6) = < VmSVmSVm >p , (45) 

here Wys, Wv2, Wc, Wmi, Wm2, Wmm, Wmcm are operators associated with the soft 
contributions, which have relatively simple form, Ey^, Ey2, E^j, Ef^j^, Ef^^j^,,, E^ are 
analytically found values of "soft" contributions to Aurec] these values and operators are 
written explicitly in the section 4. 

Remind that operators Vq + Vi , V2 , V3 are given by sums of all the tree, one-loop, 
and two-loop irreducible noncovariant graphs respectively. Hence first terms in < >p in 
(PSP - (^) are also given by sums of some noncovariant graphs (reducible ones, except 
V3). The sum of them is equal to the sum of expressions associated with all two- loop 
graphs of the recoil type. The operator corresponding to this sum is referred to as W. For 
convenience of references, the following notations are also used: sum of all second terms in 
< >p in (^) - ( ^5]) is denoted as Wo, and sum of all the soft contributions as E^ , hence 

AUrec = < W-Wo>p + E^ ; 

also, notations W', W^ are used, defined according to 

W^,W^o(p,p',7)lp=p'=.=o = / W^o(Pi,P2)- 

A sum of contributions to W', arising from all noncovariant graphs with the same 
topological structure (taking into account the difference of magnetic and Coulomb quanta), 
is equal to integrand associated with usual covariant Feynman graph (after the integration 
over zero components of the loop momenta). Hence the W' can be determined as sum of 
contributions of all covariant two-loop graphs of the recoil type, and there is no need to 
know contributions to the W' arising from separate noncovariant graphs. 

In the above the gauge used for noncovariant and covariant graphs was undermined to 
be the same. However sum of contributions of all covariant graphs to the W' is independent 
of gauge, and so in actual calculations any gauge may be choosen. This is quite natural; in 
terms of the NRQED (or the effective nonrelativistic Hamiltonian approach) < W — Wo >p 
is equal, besides some overall factor, to contribution ~ to the constant of pointlike 
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fermion-fermion interaction (if an appropriate regularization is used), which is determined 
from scattering amphtudes on the mass shelL 



4. THE "SOFT" CONTRIBUTIONS. 



In this section the general method of separation of "soft" contributions, described above, 
is applied to various terms of (0) - ( |3TD and to contributions of various graphs. 

There are some preliminary remarks concerning the calculation technique employed. 
First, every term contributing to W and Wq evidently may be replaced by its average 
over positronium polarizations (or, in other words, over directions of the full spin). This 
simplifies the calculations considerably. For contributions to W this average is found in 
program way, as 1/3 of sum over polarizations. For contributions to Wq the averaging is 
performed using the following formulas (to the r.h.s. of which terms independent of the 
full spin may be added): 

(aia) (aab) ^ ^ (ai<?2 ) (ab) , (46) 

(aid.) (a2b) (^aia2 - ^ (aiq) (a2q) ^ ^ ~ ^ ^1*^2 (ab + (aq)(bq)/q^ ) , (47) 
( (^iq) (^2q) W . . (^ik) (g2k) \ ^ W T ^ 1 (qk)^ \ 

(48) 

Second, let me describe calculation of contributions associated with retardation. These 
contributions have "ultrasoft" parts, for which virtual photon momenta ~ ma^ are essen- 
tial. These contributions are contained in < Vm >(6)5 < V2 >(6)5 < >(6) . Considering 
them, it is convenient to replace these operators by their spin-spin parts averaged over 
positronium polarizations, as described above. First nonrelativistic approximations (in the 
sense of the power expansion in momenta) for averaged contributions of graphs 2(a) and 
2(b) to V2 and are V:^^^ and VJ^*, 



yr(P,k) = / TTT^ 7T— ^1^2^ >^r(p,k,q) , (49) 
J (27r)'^ ini'^ 

^3 (P)k) = / 777^ 777^ ^ ^1^2 ^ F2 (P,k,qi,q2) , (50) 



where 



(27r)3 (27r)3 3 ' ' qfq^ 

^"^^'^'''^ = g + (/. + /.J/2 q + ih + U/2 ' ^^'^ 

yr*(p,k,qi,q2) = J^^Jf-^J^^ (g+(/,, + /,J/2) (g + (/p, + /.)/2) ^^^^ 

(notations for the momenta are shown at fig. 2); just the same contributions arise due to 
graphs obtained from graphs 2(a) and 2(b) by the time reversal. 
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Spin-spin part of Vm, averaged over positronium polarizations, is denoted below as Vg. 
First nonrelativistic approximation to l^s(p,k) is 



mO 



q+{fp + fk)/2 



(53) 



where 



Vr 



mO 



p=k=7=0 



2 an ^ ^ 
-— (Xio-a 
3 



q 



p 



It is convenient to consider < >(6)) < V2 >{6)^< V^ >(6) together and to cal- 
culate their overall contribution to Ai/^ec! this way of the calculation was employed, in 
particular, in and [||. The way of the calculation, used in the present paper is, in some 



sense, even shorter. Each of operators V^^, V{'^^, V^^^ is divided into two parts: 

\/2r(p,k) + ^'o(p,k) , 

(p,k) + V3'o(p,k) , 



P,k) 
P,k) 



-\/ret I 
^30 



where 



yret. 



P,k) 



1. 



mO 



+ 



q 



q + fp q + fk 

and V^Q^, Vr^Q* are obtained from VJ*^*, VJj^^* by replacing Y"/"^*, 



(54) 



Y{o\ where 



\rret. 
-'20 ' 



nr(p,k,q) 



p,k,qi,q2) 



1 



{q^h) {q + fk,) 
1 



+ 



2 V iq + fp)iq + fp^)iq + fp2] 

In the region where all the momenta are ~ ma 

Trret ^ -irret -t^ret 
~ ^mO 1 '-I 



+ 



{q + fp) {q + fp.) 

1 



(55) 



iq + fk)iq + fk^)iq + fk2) 



(56) 



-yret 
10 



'-2 



-'20 



If one writes V2Q, V^q as integrals like (^ , |50D , the integrands in them have powers, being 
expanded in the momenta, greater by 2 than those of Y{^^ , ¥2^^^ ; also V^q has power, 
being expanded in the momenta, greater by 2 than the power of V^^*. Due to this reason 

consist of "hard" contributions only. As for < V^r. >(q), it does 



< ^20 >(6), 



< ^30 >{6) 



contain soft contributions. However it may be easily calculated in just the way other soft 
contributions to < Kn >(6) are calculated, and is actually evaluated combined with them 
in sec. 4. 3. 

Now turn to contributions of V^o*? ^20*^*5 ^so^*- It is easy to see that their full contribution 
is equal to 0: 



< K 



ret 
mO 



> 



(6) 



+ 2<\/20 >f6) +2<V, 



(6) 



ret 
'30 



> 



(6) 







(57) 

This may be shown by a short direct calculation; however (^) is evident immediately 
due to the following reason. Consider calculation of the recoil contribution ~ a^rn? /M to 
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ground state hyperfine splitting of the "hydrogen", i.e. particle of mass m/2 moving around 
particle of large mass M. In this calculation V^^{p, k) , for instance, must be replaced by 
2m/M rjo* (P, k) ; entirely < 1/;;'=* + 2^''"* + 21^3^'^* >(6) is replaced by 
2m/ M < V^Q* + 2V{q^ + 2V^Q^ >(6), as may be easily seen. On the other hand, it is well 
known that when calculating hyperfine splitting in first order in m/M the magnetic inter- 
action may be regarded as instant, which immediately leads to 



< K^o + 2^0* + 2^30* >(6) = < KnO >(6) = . (58) 

The last point to discuss is calculation oi <Ul >(6), that is essentially the contribution 
of graphs containing Coulomb ladder with three or more loops. This contribution to Az/j.ec 
arises completely, with the accuracy required, from the region where all the momenta in 
the loops are of order ma, and can be easily evaluated analytically. There are two kinds of 
such contributions, namely, graphs having one or two magnetic quanta. The contribution 
of graphs having two magnetic quanta is 

^LMM = < V-mL Vra >(6) 

and does not depend, in the order required, on details of formalism used; for the first time 
it has been found in ||26|. Method of its evaluation, used in the present paper, coincides 



mainly with the method used in ||2^ (and involves evaluation in coordinate representation); 
the result obtained agrees with those found earlier and equals 

/ 791 \ 

EImm = ( ^ - ) ^ 0-3672 ma' . (59) 

The contribution of graphs having one magnetic photon does depend on details of 
formalism used. The value calculated in the present paper is 



El,, =< ( KnL ( Vi, + AT )) + h.c. > 



(6) 



This value was evaluated in two ways: in coordinate space (by method close to that of 
and by immediate integration over the momenta, using explicit form of L(p, k), quoted, 
for instance, in [EHl: 



-47rmV /5 1 2A - 1 r— 



where 



A 



(60) 

fpfk 



The results obtained by both methods are the same and equal to 



EIj^ = ^ ma' ^ 0.01562 ma' . (61) 
64 



Finally, addition (p9D and (|6l|) together results in 

< Ul >(6) ~ 0.3728ma^ 
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4.1. Irreducible Graphs. 

For calculation of < V3 >(6) it is necessary first to separate from V3 the term V30 which 
is an integral having the "divergency power" n = 0, and to use formula (|3^ ) where X21 
should be set equal to V30 ; the result has the form (|39| ) with 

Wv3 = V30 , Ey^ = < V30 >(6) . 

There exist only three essentially different graphs contributing to V30: these are graphs 
depicted at fig. 3 (a), 3(b) and 2(b). However contribution to V30, due to sum of graphs 
3(a) and 3(b), may be set equal to zero. The matter is that this sum in the region of 
small loop momenta reduces effectively to the graph 3(c), in which effective two-photon 
"seagull" vertex does not depend on the spins. This vertex appears as sum of subgraphs 
3(d) and 3(e) and corresponds to term aA^/(2m) (where A is the vector potential) in the 
nonrelativistic Hamiltonian of particle in magnetic field. Thus V30 arises from the only 
graph 2(b) and may be choosen to be 2^30'^* which results in Ey^ = 2 < V^q^ >(6) ; in fact 
there is no need to evaluate this value, as mentioned above. 

Similarly < V2 >(6) is calculated using (|38|) with Xn = V20 where V20 is part of V2, 
having n = . Noting that < Uv2 >(6)=< Uv2 >p , as is easily seen, one finds ( POD where 

Wv2 = {V20 + {V2 - 1^20)0 ) SoVo , = < ^20 >(6) • 

Consider contributions to ¥20- Graphs having two Coulomb quanta do not contribute to 
¥20- Among graphs having two magnetic quanta there are those with n = —1 (graphs 
depicted in fig. 3(f), (g)), and several graphs with n = 0. However their contributions to 
V20 cancel each other, just as described above for V30; for instance, sum of graphs shown in 
fig. 3(f) and 3(g) reduces effectively in the region of small momenta to the graph 3(h), in 
which effective seagull vertex does not depend on spins. Graphs having one magnetic and 
one Coulomb quanta, contributing to V2, are divided into 4 sets belonging each to one of 
covariant graphs (A),(B), shown in fig.4, and the graphs obtained by the time reversal. It 
is convenient to consider graphs including electron-positron pairs and those without pairs 
separately. 

Contributions of graphs including pairs to V20 may be choosen to be 

V2i{p,p') = V- = ....= / 1^-1-2^ ((^^ + ^")^. - k^R^) - (62) 

(the momenta notations are shown in fig. 4, k = p' — p , and Rk = /{im? + k"^) for 
an arbitrary k ). It is easy to find their overall contribution to -Ey2? ^^^^ equals 

Ep = - ma^ In a . 
^ 3 

The only graphs without pairs, contributing to V2, are the graph shown in fig. 2(a), 
discussed above, and the similar graph obtained from 2(a) by the time reversal. Their 
contribution to V20 is choosen to be 2^20*^*. So 

E^^ = Ef, + 2< >(6) . (63) 
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The second term in ( |63D need not in fact be evaluated, as explained above. 



4.2. Graphs with one magnetic photon. 



Consider calculation of < Um >(6) ? i-e- contribution to Au^ec due to reducible graphs 
with one magnetic photon, except corresponding contribution to < Ul >{6) , considered 
above. 

As discussed above, it is convenient to calculate < K >(6) instead of < Kn >{6)) where 
Vs is spin-spin part of Vm, averaged over positronium polarizations. Let us write Vg as an 
"expansion in powers of momenta": 

Vs = v::^ + v^^ + v^2 , (64) 

where V^q is defined in (^, Vmi(p, k, 7) ~ K — at p ^ k ^ 'y <^ m , and Vmi at 
p, k,^ m equals approximately a uniform function of second power. It can be shown 
that < Vm2 >(6) can be evaluated using ( p6D with X = Vm2 ■ Actually there is no need 
to evaluate < V^q* >(6) , as mentioned above. The value < Vmi >(6) can be easily found 
in analytic form, if one uses appropriate choice for Vmi ■ It is convenient to choose 



Vml{p,k) = ^— ^ VmO 



^ (Rp + R,) - \ {p^R^ + k^R,) - (p'R^ + ^R, 



q 

(65) 



where (q = p — k) . 

In the lowest order in momenta Vm ~ Vmb where Vmb is part of the Breit Hamiltonian 



due to magnetic photon exchange. The value < Um2 >(6) is calculated using ( pS] ) with 
-^11 = Vmb{SoV(i + 5*11 Vo) + h.c. , where Su, V/i are "first terms of expansion" of S — So 
and Vc — Vq in the momenta (namely, V(i{p, k) ~ — Vq, 5*11 (p) ~ 5* — 6*0 
at p k ^ 'J <^ m , and Sn , Vd at p,k <^ m equal approximately to uniform functions 
of zero power). 

The value < Ums >(6) is calculated using (Q) with 

X21 = VmbSoVoiSoVl, + SM) + h.c. , where ^i2(p), Vl.ip, k) ^ ^ii(p), ^^(p, k) 
at p,k <^ m (but generally 5*12, VY2 7^ ^n, V/i ; the explicit way of the decomposition of 

S — So and Vc — Vq is determined by the purpose of convenience of subsequent calcula- 
tions) . 

It is easy to see that in Xu, X21 defined in such manner it is enough to take into 
account only spin- independent parts of V(i, VI2 (denoted as Vii, V12 respectively) and spin- 
spin part of Vmbi that can be replaced by its average over positronium polarizations, i.e. 
VmQ- So < Um2 >{6) and <Umz >(6) can be evaluated using 

= VmoiSoVu + SiiVo) + h.c. , X21 = VmoSoVo{SoVi2 + S12V0) + h.c. . 

It is convenient to choose 
Vi2(p, k) = aTr/m^Rp = aix/m? (1 +p^/m^)~^ ; 
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Vii is defined so that Vii(p, k) = Vu(p, 0) ; under this condition an exphcit form of Vu 
does not enter final expressions. 

It is easy to find that Su, Su may be set to 

Suip) = Su{p)U=o + \j^-\j2 ^ (^=1,2); (66) 

first term in ( ]66| ) is choosen to be {S — 50)1(7=0) for Su and —Rp/A for 812- 
After a simple calculation < Um >(6) is found to have the form (|4^) with 

1 , 53 

-mo; + 

3 192 



Ei = ma^ --Ina + — + < V;;^o* >{6) , 



wmi = msv„:),SoVo+h.c.) + \^o5o(k;o +K^i) w - ((K)5o(k;o +Kni))o s,v,+h.c. 

Wm2 = KnO W(5o\42 + 5l2^) + ( Kn^^o + (Kn5(K " ^) )o ) W 

(when deriving Wm2 the term 

( (Kn — Kno)'S'llVo — {{Vm — Kno)'S'll Vo)o ) •S'oVq , 

was added to it, which causes no additional contribution of order ma^). 



4.3. Coulomb graphs. 

The Coulomb contribution due to reducible graphs to Az/^ec is equal to < Uc >(6)- 
This value was found in the style described above: < Vic >(6)? < Uc2 >{&ji < f^cs >(6) (i-e. 



separate terms of < Uc >(6)) were evaluated using (^,^,0) respectively. The calculation 
is simplified by the fact that in lowest (i.e. zeroth) order in momenta Vic contains only 
spin-independent terms, and besides that 

V^i^(p,k)|p=o = ^i^(p,k)|k=o = , 
where Vf^ is spin-spin part of Vic. 

It is easy to see that < Uc >(6) is determined by the formula (^), in which 

Wc = V.SoVlSoV, , = < VI >(6) ; (67) 

here Vl^ is first nonvanishing term of "expansion" of V^f in momenta (exactly, VYc(p, k) ^ V^ic(P) 
at p ~ -C m, and Vl^ at p,k m equals approximately a uniform function of second 
power). 

In (|67|) one can replace V-[^ by its average over positronium polarizations (denoted as 
V/c); it is convenient to set 

KiP,^) = ^-la, ^ {-U$-P'-P'^) ^r. + (P-k) + ^i?.) ; 
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using it, one easily obtains 



Er = ma^ ( \na + - 

^ 48 V 4 



4.4. Graphs with two magnetic quanta. 

Turn to calculation of < Umm2 + Ummz + Umcm >(6) , i-e. contribution of reducible 
graphs with two magnetic quanta, except corresponding contribution to < Ul >(6) , con- 
sidered above. The value < Umm2 >(6) is found using ( pHD with Xn = Vmb Sq R Vmb , 
where -R(p, p') = (27r)^(5'^(p — p') Rp. Similarly, the value < Umcm >(6) is found using 
(H) with X21 = Vmb SqVoSo R Vmb ■ It is easy to see that < Umm3 >(6) = < Umms >p ■ 
The results have form (|43|),fH) with 



WmM — Vmb So R Vmb SqVq + ( VmSVm — VmbSoRVmb )o'S'oVo , 

WmCM = Vmb SoVqSo R Vmb , 
^MM — < VmbSoRVmb >(6) ) ^MCM — < VmbSoVoSoRVmb >(6) • 

Spin-orbital part of Vmb does not contribute to -Eff^v/ ^mcm- Spin-spin and spin- 
independent contributions to Vmb are denoted as Vsb and Vw respectively; they are equal 
to 

Vsbip, P') = ^ 0^10^2 — o (^iq) (^2q) 

VibM) = -4— - + ^ ' q = p -p • 



function V{^ such that 



As the values to calculate are average values over S-state, Vib may be replaced by any 

j l^,(p,p')rf^pp' = / V^;(p,p')c?^f 



'pp ' 

It is convenient to choose 

^'(-') ( 4- H ? ^ ? ) ) ■ 

For separate parts of E^j^, Efjuj^^ one obtains 

< V;,SoRVsb + VsbSoRV;, >(6) = - 1 Ina - i ) ma^ 

< VsbS^RVsb >(6) = ( " ^ ~ ^ ^ "^"^ 



< Vl,SoVoSoRVsb + V^bSoVoSoRVl, >(6) = {~ \ + ^ 
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which results in 



/ 5 1 \ 

< VsbSoVoSoRVsb >(6) = - — Ina + — j ma^ 

t:^S ^ 3 , 23 \ 6 

Emm = ~ 8 ~ ^ J ' 

Now the calculation of the "soft" contributions is completed. Their sum is 

^ ( - - Ina + ) ma^ f - - Ina + 0.8061 ] ma^ . 

\ 6 1728 / V 6 / 

To check the formalism and the calculation technique used, they were applied to cal- 
culation of the recoil contribution ~ a^m? /M to the ground state hypcrfine splitting of 
the hydrogen. It is well known that for this particular purpose the two-particle problem 
reduces effectively to the problem of motion in external field, and the value desired may 
be easily obtained by means of coordinate-space calculations, using known solutions of the 
Dirac equation in the Coulomb field. The correction to the hyperfine splitting thus ob- 
tained is equal to 4a^m^/M. The same result was found using momenta representation 
by the method described in the present paper. To simplify corresponding calculations, 
diagram technique rules were modified, because in the relevant order one can consider the 
magnetic field as permanent, and so the magnetic photon exchange may be assumed to 
be instantaneous. All the "hard" contributions cancel each other; the soft contributions of 
the desired order arise only from graphs having one magnetic photon, and their calculation 
coincides mainly with the calculation of corresponding contributions for positronium. 



5. THE "HARD" CONTRIBUTIONS. 



For the calculation of the "hard" contribution to ^Vrec , i-e. the value < W ~ Wq >p , 
covariant graphs are used, as explained in section 3. In covariant gauge there arc only 
4 essentially different two-loop "recoil" graphs, whereas in the Coulomb gauge there are 
24 ones (regarding magnetic and Coulomb quanta as depicted by different lines); here 
term "essentially different" mean that these graphs cannot be obtained from each other by 
transposition of two particles and/or time reversal. 

The hard contribution to Aurec was evaluated in two ways. The first way is to use the 
Coulomb gauge and to evaluate contributions of the 24 graphs separately. The second one 
is to calculate < W — Wq >p as a whole using the Feynman gauge; in this case the W is 
sum of contributions of only 4 graphs. Besides that, these contributions have simpler form 
than those of the Coulomb gauge (strictly speaking, the "simpler form" refers to simpler 
form of the program calculating these contributions) . 
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The result for the "hard" contribution is — 0.424(6)ma^ for the "separate" and — 0.426(6)ma^ 
for the "united" calculation. Adding the "soft" contribution one obtains 

Ai/^ec = ma^ ( -llna + 0.382(6)) (68) 
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and 



Aurec = ma^ (^-^Ina + 0.380(6)^ (69) 

for the "separate" and the "united" calculation methods respectively. These results are in 
perfect agreement with those found in |jl|, |^ . 



The error estimation quoted in (|68| . |69| ) is arbitrary in some extent, being subject of 
handwork; the reason of this is that method used in the program for calculation of W is 
numerically unstable, and the error quoted arises from this unstability (and not from the 
error of numerical integration). The direct way of the error estimation is to compare the 
results of "separate" and "united" calculations; difference between them is indeed within 
the estimation quoted above, as well as their differences from the results of works |Q, 0, |^. 

Results for contributions of separate graphs are quoted in tables la, lb, 2 together with 
the revised results of paper , which agree with the results of the present paper. Errors are 
not given in the tables because they are not well defined; in fact, error for any contribution 
is ~ 1%. 

Notations for the graphs are following. If one makes no difference between magnetic 
and Coulomb lines, there exist 4 essentially different two-loop recoil graphs, denoted as 
1,2,3,4 at fig. 5. To denote magnetic and Coulomb lines, indices "M" and "C" are written 
in the sequence the lines are attached to the lower fermion line from left to right. Values 
quoted in the tables refer to "essentially different" graphs only, including in each quoted 
value contributions of all graphs of the similar form. 

Having been averaged over the full spin direction, W depends on three variables; con- 
venient choice for these variables are absolute values of the loop momenta, referred as 
Pi, p2, q ■ For the "united" method of the calculation it is convenient to symmetrize W 
over permutations of pi, p2i Q , and to take into account only contributions of "essentially 
different" graphs, taking some of them multiplied by 2 or 4 because of transposition of the 
two particles and/or time reversal. 

Let us describe briefly the checking of the calculation of the "hard" contributions. 
The first method of the checking is to compare < W — Wo >p found in the two ways, the 
"separate" and the "united" (and intermediate comparison of W{pi,p2, q) obtained in the 
two ways, at several points in pi,p2, q)- 

Second, contributions to Wq were found by a relatively short analytic calculation (except 
contributions associated with sum of graphs 1-CMM and 2-CMM), whereas contributions 
to W are results of the numerical computing. An error in W or Wq leads usually to 
divergency of the numerical integration; clearly search for such errors is more convenient 
for the "separate" calculation. 

Third, the results of the "separate" calculation are compared with the results of in 
which Az/j-ec was calculated using Bethe-Salpeter formalism and the Coulomb gauge. The 
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correction to the energy can be written as an expansion in irreducible covariant graphs 
(except the tree graph with Coulomb photon). Numerical value of any term of this ex- 
pansion is clearly independent of formalism used, e.g. the formalism of the present paper 
or that of 0. So most of contributions found in the present paper should be equal to 
corresponding contributions found in . The other terms can be divided into groups, and 
total contribution of each group should also be equal to corresponding contribution of 0. 
Besides that, to order under the consideration the contribution of graph 1-MCM should be 
equal separately to corresponding contribution of 0. 

Results for "hard" contributions of separate graphs are given in tables la, lb. Part of 
them (namely, those quoted in table lb)) are not uniquely defined as they depend on the 
choice of Wq- On the other hand, contribution of every term of the expansion in irreducible 
graphs, mentioned above, does not depend on formalism used, e.g. on the choice of Wq. 
Some terms of this expansion are trivial, i.e. they contain only "hard" parts, and every of 
them consist of contribution of just one covariant graph; corresponding results are quoted in 
table la. The other ("non-trivial") terms of this expansion involve also "soft" contributions, 
and their "hard" parts may consist of contributions of several covariant graphs; their values 
are given in table 2. Tables la,2 also contain corresponding results of 0]. All the results 
of are in good agreement with those of the present work. 

Contributions of graphs 3-MMC and 4-MMC were computed together, because these 
graphs contain noncovariant graphs shown in fig. 3(a) and 3(b), contributions of which to 
V30 cancel each other, which simplified the calculation. Contributions of graphs 1-CMM 
and 2-CMM were also evaluated together due to the similar reason. 

In table 2 the value E^^^q, for instance, denotes the "hard" part of the contribution of 
graph 1-CCC, i.e. one of values quoted in table lb. The value Em stands for 

EfcCM + ^2CMC + ^lAICC + ^ICMC + ^LM + + ) 

and the E'^ denotes 

(ctO+tcO)+ctcO+(cte+tcx) + (cctO+tccO) + (ctcy+cct?/)+tec2+AEA/p((5Js:o-^+^-'5^o)+A^d^' 
(using notations of paper 0). 
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Table la. 

"Trivial" (i.e. containing only the "hard" part) contributions to Ai/, 
(in units ma^ ); the results of are also quoted. 



Graph 


contribution to Az/,.ec 


(notations of 0) 


contribution to Az/,.ec 
(obtained in |Q) 




-U.U0o9 


cccx 


-0.0039 




(J.UU42 


ctcx 


r\ r\r\ AO 

0.0043 


4-MuL 


-0.0486 


cctx + tccx 


-0.0489 


4-MCM 


-0.0230 


tctx 


-0.0230 


1-MMM 


0.0o94 


ttto 


0.0694 


O A /TA /TA /r 

z-MMM 


-0.0011 


ttty 


-0.0011 


O A /TA /TA /r 

o-MMM 


-0.0012 


tttz 


-0.0011 


A A /TA /TA /T 


r\ f\r\ Ac\ 

0.0042 


tttx 


n nn /i 1 

0.0041 




0.00d4 


cccz 


n nn^ o 

0.0063 




0.02DO 


ctcz+cctz 


0.0263 


3-CMM 


0.0530 


ttcz+tctz 


0.0534 


2-CCC 


-0.0184 


cccy+ccx 


-0.0186 


2-MCC 


-0.0661 


tccy 


-0.0681 


2-MCM 


0.0552 


tcty+ttcy 


0.0558 
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Table lb. 



"Hard" parts of "non-trivial" contributions to Aurec 


(in units ma^ ). 


Graph 


1-CCC 


-0.0094 


1-MCM 


-0.0745 


1-CMM + 2-CMM 


-0.5112 


4-MMC + 3-MMC 


-0.0209 


3-CCM 


-0.0104 


2-CMC 


-0.0245 


1-MCC 


0.1795 


1-CMC 


-0.0092 



Table 2. 

"Non-trivial" contributions to Az/^ec (in units ma^ ), 
compared with the results of 0. 



value quoted 


value quoted 
(notations of 0) 


coefficient 
at In a 


constant 


constant 
(obtained in [0]) 


^4MMC+3MMC 


(cttx-t-ttcx) +cttz 





-0.0209 


-0.0209 


^ICCC + 


ccO+cccO 


-1/48 


-0.0146 


-0.0148 


E\MCM + E'mcM 


tctO 


-5/48 


0.3176 


-0.3138 


EiCMM+2CMM + Emm 


cttO+ttcO+ctty + 
+ ttO+ttx 


-3/8 


-0.751 


-0.749 


Em 


E'm 


1/3 


0.427 


0.423 



6. SUMMARY AND DISCUSSION. 



In the present paper the "recoil" contribution of order ma^ to the hyperfine splitting 
of positronium ground state has been found. The calculation was performed using nonco- 
variant formulation of QED perturbation theory. The result equals 



Aurec = ma^ ( -ilna + 0.381(6) 
6 
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and is in perfect agreement with the results of |jT|, ||, ^ . Let us note that methods used in ^ 
and in the present work are similar, whereas that of is completely different. So Aurec 
may be considered to be firmly established. However agreement with the experimental 
data is poor. Combining nonlogarithmic part of the recoil contribution with all other 
contributions of the relevant order (they also may be considered to be reliable), one obtains 
the complete nonlogarithmic contribution ~ ma^ to Au: 

ma^ ( -0.3928 ) = - 7.33 MHz (70) 

((|70D was obtained using for Aurec the result of 0, as the most accurate). The complete 
theoretical result up to the order ma^ is sum of ( [7D| ) and (^, and equals 

Auth = 203 392.96 MHz . 

It differs from the experimental result (|^) by 5 standard deviations. 

The leading term of the next order in a, i.e. of order ~ ma'' In^ a , was found in [28 
and equals 

7 

ma^ln^a = -0.92MHz . 

871 

Taking this correction into account reduces the difference between experimental and the- 
oretical results to 4 standard deviations. One cannot reject possibility that this difference 
may be explained by contributions ~ ma^ In a. 



ma 



7 
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Fig. 5 
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